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We study s-wave superconductivity in the two-dimensional square lattice attractive Hubbard 
Hamiltonian for various inhomogeneous patterns of interacting sites. Using the Bogoliubov-de 
Gennes (BdG) mean field approximation, we obtain the phase diagram for inhomogeneous patterns 
in which the on-site attractive interaction Ui between the electrons takes on two values, Ui = Q 
and —U/{1 — f) (with / the concentration of non-interacting sites) as a function of average electron 
occupation per site n, and study the evolution of the phase diagram as / varies. In certain regions 
of the phase diagram, inhomogeneity results in a larger zero temperature average pairing amplitude 
(order parameter) and also a higher superconducting (SC) critical temperature Tc, relative to a 
uniform system with the same mean interaction strength {Ui = —U on all sites). These effects are 
observed for stripes, checkerboard, and even random patterns of the attractive centers, suggesting 
that the pattern of inhomogeneity is unimportant. The phase diagrams also include regions where 
superconductivity is obliterated due to the formation of various charge ordered phases. The en- 
hancement of Tc due to inhomogeneity is robust as long as the electron doping per site n is less 
than twice the fraction of interacting sites [2(1 — /)] regardless of the pattern. We also show that 
for certain inhomogeneous patterns, when n = 2(1 — /), increasing temperature can work against 
the stability of existing charge ordered phases for large / and as a result, enhance Tc. 

PACS numbers: 



I. MOTIVATION 

Fascination with inhomogeneous superconducting 
(SC) phases extends back several decades, with many 
conferences and monographs having been devoted to the 
subject jii^ii^ For conventional superconductors, the inho- 
mogeneities were extrinsic, arising from a granular na- 
ture of samples or due to the deliberate synthesis of dis- 
ordered built materials or films. The high temperature 
superconductors (HTS) introduced new aspects into this 
area of study. First, inhomogeneous states (normal and 
SC) seem to be intrinsic to HTS, at least in the under- 
doped regime, similar to quenched disorder in the metal- 
insulator transition in two dimensionsi^ Secondly, the 
inhomogeneity occurs on a smaller length scale of just a 
few lattice constants as evidenced by scanning tunneling 
spectroscopy (STS) at the nanoscalei^i^ii^ 

In addition, the strong electronic correlation in HTS 
cuprates plays a major role in the elucidation of the inho- 
mogeneous SC state; indeed the inhomogeneity is widely 
believed to arise from the strong intra-atomic interactions 
that tend to frustrate bandlike conduction, to induce lo- 
cal magnetic moments, and to drive charge and spin order 
on a few lattice constant scale. The study of the doped 
Mott insulating phase has been one of the most active 
areas of theoretical study of HTS, one that has proven 
to be unexpectedly complex and rich and which is still 
under intense exploration and debate. Surprisingly, holes 
doped into the high temperature superconductors tend to 
arrange themselves non-uniformly in the Cu02 planes in 



the form of stripes, checkerboard or perhaps even more 
exotic structures . "'^^d^d^d^, 14 ]y[oreover, spatially varying 
density and spin structures have also been observed in the 
physics of the manganite o-'^^d^d'^d^ and cobaltitesji^i^St 

A variety of physically relevant models such as 
the repulsive Hubbard and t-J Hamiltonians have 
been extensively studied to understand the inter- 
play between spatial inhomogeneity, magnetism, and 
^upcrconducti'.'ity .^^-^^'^'^'^'^-^^-^^-^'''^^-^^''^"-^^-'^^-'^'^''^'^-^^-'^^-^''-^^-'^^-^" 
In the repulsive Hubbard and t-J Hamiltonians in partic- 
ular, inhomogeneity has been introduced either through 
the hopping amplitude t or magnetic coupling J or the 
local energy on the lattice sites. For the 2D square 
lattice these two models are known to display antiferro- 
magnetism at half-filling, and, although it is less certain, 
perhaps also d-wave superconductivity when doped. 
There is considerable evidence that they also might 
possess inhomogeneous stripe or checkerboard ground 
states;^i2^iS2i2^i2ai2ii2^i2i. Phenomenological d-wave 
BCS Hamiltonians with spatially inhomogeneous pairing 
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have also 



been employed mostly to reproduce the local density of 
states (LDOS) results obtained from scanning tunneling 
microscopy (STM)4^ In addition, there have been 
theoretical studies of the SC quantum phase fluctuations 
using the QED3 effective theory of the HTS in relation 
with the inhomogeneous pattern formation in cuprates 
from the STM experimental results.*^'— 

While DMRG treatments'^'* provide detailed informa- 
tion on the real space charge, spin, and pairing orders, 
the precise nature of the interplay, and whether the dif- 
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ferent orders compete or cooperate, remains unclear. In 
addition, the enhancement of the superconducting tran- 
sition temperature Tc by local inhomogeneity has been 
demonstrated by Martin et al. in Ref. [48^ and also in 
Ref. i49j]- Recently, the XY model Hamiltonian with 
certain types of inhomogeneous patterns for the coupling 
constant between spins sitting on two nearest neighbor- 
ing sites has also been shown to enhance Tc by Loh et al. 
in Ref. 

Many of the basic characteristics of this short-range- 
scale inhomogeneous superconducting state can be ad- 
dressed with a more tractable model, one which is well 
understood in the homogeneous limit. This model is the 
attractive Hubbard model, which has been applied pre- 
viously to address some aspects of the impact of inho- 
mogeneity. Recently old suggestions based on "negative 
U" superconductivity have been revrved^^ which may 
provide additional applications for the results we present 
in this paper. Tl-doped PbTe achieves a critical tem- 
perature up to 1.5 K, and more extensive heat capacity 
and transport data^i^ have led to an analysis in terms 
of a "charge Kondo effect" that could be linked to the 
observed superconductivity.^"* This system intrinsically 
involves both negative U centers and inhomogeneity. 

This article extends previous work^ to a more general 
range of non-interacting site concentration / values. We 
show the presence of different conduction phases both 
in the phase diagram at T = and in the density of 
states (DOS). We also show the local occupation and 
SC order parameter for electrons on different lattice sites 
as the concentration / varies for different inhomogeneity 
patterns. Finally, the Tc enhancement conditions are also 
extracted based on the relationship between the average 
doping of electrons n on the lattice and inhomogeneity 
concentration /. 

The paper is organized as follows: in the next section 
we introduce our model and describe the method we have 
employed. In Sec. Ill we present and discuss the phase 
diagram at zero temperature. Sec. IV contains our finite 
temperature results, and Sec. V summarizes our findings. 



II. MODEL AND METHODOLOGY 

This article focuses on a general question: "Under 
what conditions is it more favorable to have an inhomo- 
geneous pairing attraction, compared to the same aver- 
age strength spread homogeneously over the lattice?" By 
'conditions' we refer to the average occupation number 
of electrons per site n, the average attraction strength 
per lattice site U , which remains constant in any com- 
parison between systems with and without inhomogene- 
ity, and the inhomogeneity concentration /. We address 
this question by comparing the average zero temperature 
pairing amplitude A over the entire lattice and the SC 
transition temperature for a system in the presence 
and absence of inhomogeneity. 

For the cuprate superconductors, for example, such a 



question is complicated by the presence of other types of 
order such as antiferromagnetism, exotic spin-gap phases, 
and by the nontrivial d-wave symmetry of the SC or- 
der parameter. For these systems and phenomena, mod- 
els like the repulsive Hubbard or t-J Hamiltonians are 

^,„„j,^|-^^^ 21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36.37.38 jv^gy- 

ertheless, it is yet beneficial to study the problem first by 
employing a more simple and phenomenological model. 
Here we will present a solution of the inhomogeneous 
Bogoliubov-de Gennes (BdG) equations for the attrac- 
tive Hubbard Hamiltonian, 

<ij>,a 

- tJ''Y4aCta-Y\U.,\ni^nii, (1) 

i(j i 

with t the hopping amplitude, fi the chemical poten- 
tial and Ui the local attractive interaction between the 
fermions of opposite spins residing on the same lattice site 
i. We will study inhomogeneous patterns in the interac- 
tion Ui. The origin of the attraction in this model can 
result from, for example, integrating out a local phonon 
modci^ The two-dimensional uniform attractive Hub- 
bard model is known to yield degenerate superconduc- 
tivity and charge density wave (CDW) long range or- 
der at half-filling and zero temperature fSSiSi^ However, 
away from half-filling, the CDW pairing symmetry is bro- 
ken and superconductivity is more favorable, and the SC 
phase transition is at finite temperature. 

The BdG mean field decomposition bilinearizes the 
Hamiltonian by replacing the local pairing amplitude 
and local density by their average values. A; = (ci^Cii) 

and (riiij) ~ (^c^^^Cia) and yields the quadratic effective 
Hamiltonian 

<ij>.(T 1(7 

i 

where jli = /i + |C/i|(ni)/2 includes a site-dependent 
Hartree shift with {rii) — '}2,^{nic). All energies will be 
referenced to t = 1. 

We adopt the criterion of comparing the tendency for 
superconductivity in the homogeneous system with the 
same attraction —[/on all lattice sites, with cases when 
sites with attraction are mixed with sites where the at- 
traction is absent, i.e., Ui = ,48,59,60 Specifically, we 
have studied systems in which sites with attractive inter- 
action are randomly distributed^ or arranged in checker- 
board and stripe patterns. The last two regular patterns 
have been purposely chosen due to their relevance to the 
experimentally observed pattern formation in the HTS 
cuprates. 

Fig. [U panel (a), presents the patterns for the in- 
teracting lattice sites with four different values for the 
fraction of non- interacting sites /. The uniform pattern 
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f = 0.0 
(uniform) 



f = 0.25 
(checkerboard) 



f=0.5 
(a) (checkerboard) 



f=0.75 
(checkerboai'd) 
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f= 0.25 and 0.75 
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IUi = (4/3)U 
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f = 0.25 



U. = 2U 



U:=0 

f = 0.5 



Ui = 4U 

U2=0^ 

f = 0.75 



FIG. 1: (Color online) Panel (a): Regular patterns for the in- 
teracting sites in the attractive Hubbard model at diflFerent in- 
homogeneity concentration values. Stripes and checkerboard 
have been particularly selected because of their relevance to 
the experimental observations in cuprates. Panel (b): Color 
coding and numbering different types of sites for the checker- 
board and stripes blocks as presented in panel (a) (two colors 
for / = 0.5 and four for / = 0.25 and 0.75). Sites carrying 
identical color code and number are equivalent by the sym- 
metry in the lattice geometry. Panel (c): Color coding of the 
lattice sites for the random inhomogeneous pattern at differ- 
ent / values. Regions of interacting (non-interacting) sites 
are coded black (red or dark gray in the grayscale version) 
with the appropriate weight of 1 — / (/). 



corresponds to / = with interaction Ui — U — 
on all lattice sites. Checkerboard, stripes and random 
patterns with / = 0.25 include 1 — / = 75% interact- 
ing sites with Ui — |t/ = — and / = 25% non- 
interacting sites with Ui ~ 0. For / = 0.5, half of the sites 
are interacting with Ui = —2U and half non-interacting 
with Ui = 0. / = 0.75 corresponds to 25% interacting 



sites with Ui = —4/7 and 75% non-interacting sites with 
Ui = once again averaging to C7 = —U per site. 

For the random pattern we have averaged over typ- 
ically 20 different disorder realizations. One may note 
that regardless of the pattern and the impurity fraction, 
the average attraction per site, i.e., U = —U remains 
constant. We adopt this criterion throughout the rest of 
this article for comparison between uniform and inhomo- 
geneous lattices. This criterion is believed to be most 
appropriate for exhibiting the effect of inhomogeneity in 
particular in the systems having the same strength on 
average for forming superconducting Cooper pairs. 

Panel (b) depicts the color coding and numbering of 
the lattice sites for the checkerboard and stripes based 
on their value of /. Due to the regular geometry and pe- 
riodicity of their inhomogeneous patterns, lattices with 
the checkerboard and stripe patterns can be classified 
into blocks including two and four different types of sites 
at / = 0.5 and / ~ 0.25 or 0.75 respectively as illustrated 
by different colors and numbers in panel (b) of Fig. [T] 
Sites carrying the same color or number are equivalent by 
the symmetry in the pattern geometry. For the random 
pattern, due to the lack of both regularity and periodic- 
ity, we average over all the interacting and all the non- 
interacting sites separately (black and red (Red shows 
dark gray in the grayscale version) regions in panel (c) of 
FigU respectively) before the configurational averaging 
over all different impurity patterns is performed. 

We self-consistently diagonalize the BdG mean-field 
Hamiltonian in Eq. [2] by assigning initial values to the 
local occupation number rii and order parameter and 
solving again for these parameters after diagonalizing the 
Hamiltonian, until convergence is achieved at a desired 
tolerance. For the checkerboard and striped patterns, 
sites with the same color have the same density, due to 
symmetry, and do not need to be averaged. For the ran- 
dom pattern we calculate the average occupation number 
?^coior and order parameter Acoior per site at the different 
types of sites by averaging 



^color 
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For the checkerboard and striped patterns with / = 
0.5 and for the random pattern (all /) we simply have 
"attrac = "black and rifree = "rod (Red shows dark gray in 
the grayscale version). For the striped and checkerboard 
with / = 0.25 we have nattrac = (ngrocn + 2nrcd)/3 (Green 
shows light gray in the grayscale version) and rifreo = 
"-black, whereas for / = 0.75 we have riattrac = "black and 
"free — ("green + 2nred)/3. The Same Combinations hold 
for Aattrac and Afree- The average occupation number n 
and order parameter Ajnhom per site are defined 



Ai 



= (1 - /)"attrac ' 
(1 ~ /)Aattrac " 



/"free 
/Afree 



(4) 
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The chemical potential in Eq. [2] is self-consistently 
adjusted after every iteration in order to arrive at a de- 
sired total average occupation per site n for the entire 
lattice. For the regular patterns, i.e., uniform, checker- 
board and stripes, due to their periodicity, by Fourier 
transforming the Hamiltonian into momentum space, we 
significantly reduce the numerical cost of the calculations 
and at the same time can increase the size of the lattice 
close to the thermodynamic limit to avoid finite size arti- 
facts in the results (up to 1500 x 1500 lattice sites in our 
calculations). For the random pattern, however, such 
a simplification is not possible due to the lack of peri- 
odicity. Hence, we are limited to the finite size lattices 
of up to 24 X 24 sites. As a result, especially at small 
values for the average on-site interaction magnitude \U\, 
finite size effects are to be cautiously monitored. Our 
calculations also include the density of states (DOS) for 
the entire lattice. We study simultaneous variations of 
the average on-site interaction magnitude \U\, occupa- 
tion number n and also temperature and their effects in 
the average order parameter Ainhom and DOS. Our goal 
is to obtain the phase diagram for the effect of inhomo- 
geneity in superconductivity and discuss the conditions 
under which inhomogeneity can result in enhancements 
in the average superconductivity order parameter or SC 
phase transition temperature Tc- 

It is further realized that our conventional mean-field 
approach does not capture the Kosterlitz-Thouless na- 
ture of the phase transition in two dimensions. Never- 
theless, this weakness can be repaired^ ^ upon regarding 
the local pairing amplitudes as complex variables and 
performing a finite temperature Monte Carlo integration 
over the associated amplitude and phase degrees of free- 
dom. Unhke BCS, this Monte Carlo mean field (MCMF) 
approach allows identification of the weak and strong 
coupling regimes via the phase correlation function. In an 
earlier work*^ this Monte Carlo technique was employed 
as an independent examination for the validity of our re- 
sults and the agreement between the two techniques was 
clearly confirmed. 



III. PHASE DIAGRAM AT T = 

Fig. [2] presents the phase diagram for the average 
interaction magnitude \U\ and electron doping n per 
site at T = for three different inhomogeneous pat- 
terns of checkerboard, stripes, and random, and for / 
values of 0.25, 0.5 and 0.75. We show isocontours of 
r = Ainhom/ Auniform, 1-6., the ratio of the average inho- 
mogeneous pattern order parameter over its uniform pat- 
tern counterpart. Thus, whenever r > 1, inhomogeneity 
leads to a larger average order parameter at T = com- 
pared to a homogeneous system and therefore is more 
favorable for superconductivity over a uniform pattern 
of the interacting sites. We also adopt the lower limit of 
= t since for smaller values of \U\, r will be the ratio 
of two very small numbers and is subject to numerical 



uncertainty. The first row of Fig. [2] (panels (a)-(c)) cor- 
responds to the concentration value oi f — 0.25 for the 
non-interacting sites. At first glance, one can observe 
that regardless of the geometry for the inhomogeneity, 
above n — 1.5 = 2(1 — /), inhomogeneity gradually (or 
abruptly for the checkerboard pattern in panel (a)) re- 
sults in the obliteration of superconductivity consistent 
with the findings of Litak et a^. We can understand how 
this obliteration takes place if we examine the behavior 
of the system in strong coupling. When we start with an 
empty system and add electrons they are placed on the 
attractive sites first due to the strong attractive interac- 
tions. It is useful to define n* =2(1—/), which for a given 
/ corresponds to the density for which all attractive sites 
are doubly occupied and all free ones are empty. As we 
will see below, this density corresponds to an insulating 
charge ordered state. Above this density, superconduc- 
tivity is reduced because the pairs cannot move within 
the attractive sublattice, since it is completely filled. 

For the checkerboard pattern in panel (a), there are 
two insulating regions within the phase diagram at n = 1 
and 1.5 (hatched orange lines) both corresponding to the 
formation of charge ordered phases of electrons in the 
interacting sites. No superconductivity was observed for 
n ~ I and 1.5 down to the lower limit of \U\ — t. Beyond 
n = 1.5, the system becomes metallic. For stripes as 
shown in panel (b), similar features as in panel (a) are 
observed. One exception is the lack of the charge ordered 
insulating phase at n = 1. This can be the result of 
further overlap between the Cooper pairs since for the 
stripes, nearest neighboring sites are interacting in one 
dimension. The charge ordered insulating phase at n — 
1.5 also forms for rather higher \U\ values compared to 
the checkerboard pattern. The random pattern in panel 
(c) also shows features similar to the stripes. 

The second row (panels (d)-(f)) corresponds to f — 0.5 
with rather similar features to the first row. As an- 
ticipated, superconductivity gradually goes away above 
n = 1 = n* for all the patterns. For the checkerboard in 
panel (d), and n > 1, superconductivity strictly goes 
away and the system turns metallic. For the striped 
and random patterns however, superconductivity persists 
slightly above n — 1 until it is totally obliterated. At 
n = 1, all three inhomogeneous patterns exhibit a charge 
ordered insulating phase for large enough values of |C7| 
(or the smallest value of |C/| for the checkerboard). Nev- 
ertheless, it can be readily seen that for / = 0.5 compared 
to / = 0.25, the enhancement of the average order pa- 
rameter due to inhomogeneity is considerably larger as 
the ratio r increases for small \U\ values. 

The third row (panels (g)-(i)) for / — 0.75 shows ratios 
as large as r = 15 for small values of |C7| and n values not 
much larger than 0.5 = n*. For / — 0.75 also, supercon- 
ductivity gradually dies away when n > 0.5 ~ n* and a 
charge ordered phase sets in for large enough \U\ values 
at n = 0.5 = n* analogous to / = 0.25 and 0.5. The only 
difference is a slight remnance of superconductivity for 
the checkerboard pattern at n > 0.5. Thus, by further 
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diluting the interacting sites in the lattice and keeping 
the attractive pairing energy constant at the same time, 
superconductivity is driven towards smaller n values; on 
the other hand, the average order parameter becomes 
significantly more enhanced due to inhomogeneity. Gen- 
erally, regardless of the pattern, for large enough \U\ val- 
ues, inhomogeneity weakens superconductivity for every 
value of n due to the localization and compression of 
the Cooper pairs in the interacting sites. For n < n*, 
Ainhom increases as a function of |t/| and saturates for 
large \U\ values. For n > n* , Ainhom reaches a max- 
imum as a function of |C/| and eventually vanishes for 
large enough \U\ values. However, Aunitorm is a mono- 
tonically increasing function oi\U\ and is symmetric with 
respect to n = 1. Therefore, for sufficiently large \U\, the 
ratio r = Ai„hom/ Auniform bccomes less than one while 
n < n* and eventually zero when n > n* as illustrated 
in Fig. H 

The obliteration of superconductivity is associated 
with the vanishing of average order parameter Ainhom = 
0. Whether a non-SC state is a metal or insulator is de- 
termined by the DOS results for that state. In Fig. [3l 
panel (a) presents the DOS results for the checkerboard 
pattern at \U\ — 6t (the largest in our calculations for 
/ = 0.25) and four different values for the average elec- 
tron doping n. There is a gap in the DOS around the 
Fermi energy at a; = at n = 0.4. This gap corresponds 
to a SC state as for n = 0.4, Ainhom ^ 0. The gaps at 
n = 1 and n = 1.5 both correspond to insulating states 
as for both these n values Ainhom — 0. At n = 1.8, 
we also find Ainhom — according to Fig. [21 panel (a). 
However, the DOS dX n ~ 1.8 has a finite value around 
the Fermi energy at w = as shown in panel (a) of Fig. 
[3l indicating a metallic state. 

In panels (b) and (c) in Fig. [3] for striped and random 
patterns respectively, gaps at n = 1 do not correspond to 
insulating phases as opposed to panel (a) whereas they do 
correspond to insulators for n = 1.5. In the second row 
of Fig. [3] (panels (d)-(f)), for all the patterns at / = 0.5 
and \U\ = 5t (the largest in our calculations for / = 0.5), 
the system is superconducting for n < 1, insulating for 
n — 1 and metallic for n > 1 in conjunction with the 
results of the second row in Fig. [2l The third row of Fig. 
[3] (panels (g)-(i)), for / = 0.75 and |C7| = At (the largest 
of our calculations for / ~ 0.75), confirms the results 
presented in Fig. [21 (panels (g)-(i) respectively), namely, 
superconductivity for n < n* , insulator at n = 0.5 = n* 
and large enough \U\ and metal for n = 0.65 > n* for all 
inhomogeneous patterns. 

The insulating state for the phase diagram in Fig. [H 
at all values of /, n and all inhomogeneous patterns 
is always associated with the formation of a charge or- 
dered state. In Fig. [H panel (a) for the checkerboard at 
/ — 0.25, for the insulating phase at n = 1, as \U\ in- 
creases, electrons form pairs in the interacting sites with 
higher symmetry [sites 2 and 4 in red (dark gray in the 
graycsale version)], leaving the non-interacting and lower 
symmetry interacting sites [sites 1 and 3 in black and 



green (light gray in the grayscale version), respectively] 
essentially empty. For the insulating phase in panel (a) 
at n = 1.5, the lower symmetry interacting site [site 3 in 
green (light gray in the grayscale version)] also obtains a 
pair leaving only the non-interacting site (site 1 in black) 
empty. In other words, charges rearrange themselves into 
ordered pair configurations forming an insulating phase. 
According to panel (b) in Fig. [Ijfor stripes, n = 1 does 
not correspond to an insulating phase. Panel (b) in Fig. 
[H confirms this finding as the local occupation of the in- 
teracting sites (sites 2, 3 and 4 in red (dark gray in the 
grayscale version), green (light gray in the grayscale ver- 
sion) and red respectively) for large \U\ at n = 1 does 
not approach 2 while for the non-interacting sites (sites 
1 in black) it approaches zero. 

For the random pattern in panel (c) of Fig. [J] at 
/ = 0.25, we plot nbiack and firod (Red shows dark gray 
in the grayscale version) as defined in Eq. [3l The fact 
that ribiack does not approach 2 as \U\ increases (no pair 
compression) is consistent with the lack of an insulating 
phase at n = 1 for the random pattern at / = 0.25 (Fig. 
[H panel (c)). However, at n = 1.5, both striped and ran- 
dom patterns turn insulating as a result of charge ordered 
phase formation as shown in panels (b) and (c) in Fig. [4l 
where their occupation numbers on the interacting sites 
all approach 2 as \U\ increases. Charge ordered phase 
formation is more obvious in the second row (panels (d)- 
(f)) and third row (panels (g)-(i)) of Fig. [4]for / = 0.5 
and / = 0.75 respectively at large enough \U\ values. 

The proximity effect for the non-interacting sites 
neighbored by the interacting sites plays a key role in 
the magnitude of the average order parameter Ainhom in 
the inhomogeneous lattice. In Fig. [5l panel (a), for the 
checkerboard at / = 0.25 and n = 0.25, the local or- 
der parameter on all the interacting sites (2, 3, 4 in red 
(dark gray in the grayscale version), green (light gray in 
the grayscale version) and red respectively) increases as 
a function of |t7|. The non- interacting site (1 in black) 
is also superconducting due to the proximity effect of 
its neighboring sites. However, its local order parameter 
has a maximum at a critical |i7c| value beyond which it 
decreases as a result of the compression of the Cooper 
pairs in the interacting sites and therefore their weaker 
overlap around the non-interacting site. Thus, Ainhom 
on all these four sites will be larger than its uniform pat- 
tern counterpart due to this proximity effect as long as 
the non-interacting site local order parameter (Ai in this 
case) does not plummet too far down with respect to its 
maximum as a function of \U\. 

Panel (b) shows the same behavior for stripes. In panel 

(b) , there is an intersection between A2 and A3 near 
\U\ ~ 3t indicating that due to the particular symmetry 
of the stripes, sites 2 and 3 behave very closely. In panel 

(c) , we have plotted Abiack and Ared (Red shows dark 
gray in the grayscale version) as defined in Eq. [3] and it is 
clear that Aj-od eventually falls off at large |i7| values. In 
the second row of Fig. [51 (panels (d)-(f)) corresponding 
to / = 0.5 and n = 0.5, there are only two different 



6 



6 
5 
4 

l5 3 
2 
1 


6 
5H 
4 

ID 3H 
2 

H 




1 , 

p 


' 1 ' 
\ 0.75 


oi 


SC 

,' \ ^ 


'1 

(a) ; 




\ SC ; 




,' \ > 
/ , '' 


_ 

metallic ■ 


- 1.25^^ 




> 


0^5 I 

> 








-f 
< 










' 0.75 1 








-f 


^\^!25< 




n 






M 5 




checkerboard, f = 0.25 

1 1 1 1 1 1 


charge ordered insulator 

1 1 1 1 1 1 1 1 


0.2 


04 0.6 0.8 


1 

n 


1.2 1.4 1 


.6 1.8 2 



6 

5 
4 

15 3 



1.25 
1.5/- 

2, ^ 



SC \ 



metallic 



checkerboard, f = 0.5 

^ I I I I I I I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



1 

0.5 




1 1 1 1 1 1 1 1 1 


(g) 


,^SC , 


/ 
/ 


, charge ordered insulator 




0.75\ 
k \ 


\ X 










metallic 




5. 'n' 
15'~^; 


1 






checkerboard, f 

1,1, 


= 0.75 

.1.1.1.1.1. 







0.75 5 0, (b) - 


\ 


-* 




SC \ y \ 


- \ 

- 1.25~^^ ^ 


, ^\ ' \ ^ metallic . 


" 1.5 


\ * ^ 


- 1.75''~\J"~.^ 












- stripes, f = 0.25 




charge ordered insulator 





0.2 0.4 0.6 0.8 1 1.2 14 1.6 1.8 2 

n 



0.75\ 0.5 0/ (d) - 

/ charge ordered insulator " 



- P 3- 



0.75\ 0.5 


0; (e) - 


- 1 >^ 
1.25 \ 


/ charge ordered insulator - 




i 

J metallic 










- stripes, f = 0.5 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

n 



- 0.5 0/ 
5- / charge ordered insulator 



metallic 



(h) 




stripes, f = 0.75 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

n 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 l.i 

n 



ID 3 
2 
1 




1 \ 0.75 


IV ' 1 ' 1 < 1 ' 1 

05 0^ (c) - 


:\1.25 \ SC '\ 


\ f - 

1 l' 
II 

^\ ' ' 1 

\ \ I ■ metallic " 


- . I'S \ 


V II. 

^ 1 1 \ 


- 1.75^\^.^ \ 

: 


\ \ 1 1 
'\ " ^ 

\ ' \ * 
\ ^ ^ 'l 
\ \ 1 ' 


- random, f = 0.25 




1 

charge ordered insulator 

.1.1.1.1.1. 


0.2 0.4 0.6 0.8 


1 1.2 1.4 1.6 1.8 2 

n 


1 ' 1 ' 1 ' 
SC 


,,4- ' 1 ' 1 ' 1 ' 1 

»; (f) - 


' 1 0.75\ 0.5 0.25', 


J charge ordered insulator " 






- 2.J.5oy\\^ 


\ \ metalhc 










- random, f = 0.5 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


0.2 0.4 0.6 0.8 


1 1.2 1.4 1.6 1.8 2 

n 


. ' ' ' 'oi ' ' 


1 1 1 1 1 1 1 1 1 

(i) - 


SC •V 

^~ ~ - . charge ordered insulator 


: 1 0.75 0.5 ; 




- '-5 \ ^ 


metallic 






- random, f = 0.75 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

n 



FIG. 2: (Color online) Panel (a): The contour plot phase diagram for the checkerboard pattern at f — 0.25 and T — 
0. The horizontal axis presents the average occupation of electrons per site n and the vertical axis refers to the average 
interaction magnitude between two electrons per site \U\. Lines with numbers next to them correspond to different ratios of 
' — ^checlcerboard /Auniform- Solid hues at r = 1 determine the enhancement boundary. Dashes along r = lines indicate charge 
ordered insulating phase behavior. Dotted lines at \U\ = t are lower limits for the interaction as for too small \U\ values, r is 
an ill defined quantity. Panel (b): The same results for stripes with r — Astripcs /Auniform at / = 0.25 and T — 0. r — line 
for stripes is diverted towards larger n values at smaller \U\ and does not run down to arbitrarily small \U\ values at n = 1.5. 
Dashes along r = line for stripes appear only at n = 1.5 and beyond that, r — corresponds to a metallic phase. Panel 
(c) presents results for the random pattern. Similar to the stripes, for the random pattern, r = line does not run down to 
arbitrarily small \U\ values at n = 1.5 either. Panels (d), (e) and (f) correspond to / = 0.5. Note that the charge ordered 
phases for the striped and random patterns at / = 0.5 again occur only at n = 1 portion of the r = line and beyond that, 
r = yields a metal. Similarly in panels (g), (h) and (i) corresponding to / = 0.75, all three different inhomogeneous patterns 
have a r = line portion above n = 0.5 at which the systems become metallic. Panel (i) also lacks the r = 15 contour due to 
the finite size effect uncertainties at small \U\ values. 



sites for each pattern and the lattice has a more dilute 
interacting pattern. As a result, compared to / = 0.25 
results, Ainhom at / = 0.5 tends to saturate for large 
\U\ values for all the patterns. In the third row of Fig. 
[5] (panels (g)-(i)) for / = 0.75, the lattice is even more 



dilute in terms of interacting energy. Therefore, Ainhom 
shows even faster saturation at smaller \U\ values. 

As mentioned earlier in this section, for n < n* , Ainhom 
increases as a function of \U\ and saturates for large 
enough \U\ values while Auniform is a consistently in- 
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FIG. 3: (Color online) Panel (a): Density of state (DOS) for the checkerboard pattern at / = 0.25 and = 6t (the largest 
in our calculations for / = 0.25) for different values of the average electron occupation per site n. Panel (b) and panel (c): 
The same results as in panel (a) for striped and random patterns respectively. Panel (d): Results of panel (a) at / = 0.5 and 
\U\ = 5t (the largest in our calculations for / = 0.5). Panel (e) and (f): The same results as in panels (b) and (c) for striped 
and random patterns respectively. Panels (g), (h) and (i) correspond to f — 0.75 and \U\ — it (the largest in our calculations 
for / = 0.75) for the checkerboard, striped and random patterns respectively. The particular selection of colors is for better 
visibility in both the color and grayscale versions. 



creasing function of and is symmetric with respect 
to n = 1. Therefore, for n < n* as illustrated in Fig. [1] 
for r = Ainhom/Auniform at a given n 

U2>Ui^ r{U2) < r{Ui) . (5) 

Now since 

Ainhom(C^2) > A 

inhom 

Ai„hom(r(L/2)) > Ainhom(r-(C/l)) 

for [n < n*) , (6) 



and as a result 

Ai„hom(?' = 1) > Ai„hom(»' > 1)) 

for {n<n*). (7) 

Therefore, knowing that r — I yields the largest mag- 
nitude of Ainhom that is still enhanced compared to 
Auniform whcn Ti < Ti* , the optimum effect due to in- 
homogeneity corresponds to a value of n along the r — 1 
contour in Fig. [1] for which Ainhom is maximized. Fig. 
ini demonstrates the variation of Ainhom = 1) as a func- 
tion of n for different / values. In panel (a) corresponding 
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FIG. 4: (Color online) Panel (a): The evolution of the local electron occupation number Ui on different lattice sites (As color 
coded inside the blocks in Fig[T] panels (b) and (c)) for the checkerboard pattern a.t f — 0.25 and n — 1 and n = 1.5 (referring 
to charge ordered phases in Fig[2] panel (a)) as a function of \U\. Panels (b) and (c): The same results as in panel (a) for the 
striped and random patterns respectively. For the random pattern, data are taken by averaging the occupation number over 
the interacting and non interacting sites. Panels (d), (e) and (f): The same results at / — 0.5 and n = 1 (charge ordered phases 
in Fig[2l panels (d), (e) and (f)). Also panels (g), (h) and (i) at f = 0.75 and n — 0.5 (charge ordered phases in Fig[2] panels 
(g), (h) and (i)). 



to f — 0.25, all three different inhomogeneous patterns 
yield the maximum Ainhom(?' = 1) within the range of 
n = 0.5 to 0.75. In panel (b) for / = 0.5, the maxima 
are closer to n = 0.5 while in panel (c) for / — 0.75 they 
are around n = 0.25. These results indicate that ap- 
parently the optimum value for the doping of electrons 
in these inhomogeneous systems is Uopt ^ I — f = n* /2. 
In strong coupling this density corresponds to leaving the 
free sites empty and singly occupying the attractive ones. 
By comparing this behavior with the uniform system for 
which, due to particle-hole symmetry, Uopt — 1 we can 
understand why Uopt ~ n*/2. 



Also, one observes in Fig. [5] that by further diluting 
the interacting sites in a lattice, the magnitude of Ainhom 
at Hopt decreases. 

We conclude in this section that by further diluting the 
density of interacting sites in a lattice while maintaining 
the average pairing energy per site constant at T = 0, 
the average order parameter Ainhom may enhance. This 
enhancement results from the proximity effect in the non- 
interacting sites due to their interacting neighbors lead- 
ing to a larger average order parameter compared to the 
uniform lattice and in many respects is independent of 
the particular inhomogeneous pattern. Superconductiv- 
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the blocks in Fig[l] panels (b) and (c)) for the checkerboard pattern at / = 0.25 and n = 0.5 as a function of \U\. The curve 
symbolized with stars refers to the average order parameter Ainhom- Panels (b) and (c): The same results as in panel (a) for 
the striped and random patterns respectively. One notices that curves for A2 and A3 in panel (b) intersect around \U\ ~ 3t. 
Panels (d), (e) and (f) present the same results at / = 0.5 and n — 0.5 for the checkerboard, striped and random patterns 
respectively. Also panels (g), (h) and (i) at / = 0.75 and n — 0.25. 



ity in an inhomogencous lattice of interacting sites grad- 
ually vanishes above n = n* and for large enough \U\ 
values at n — n* we have an insulating phase as a re- 
sult of a charge ordered phase formation. For larger / 
values, the enhancement ratio r = Ai„hom/Auniform will 
be larger for small values of |C/| and n. However, the en- 
hancement of Ainiiom occurs at smaller values of n. We 
also find an optimum value of Uopt ~ 1 — / — n* /2 for 
the largest enhanced Ainhom for a system in the presence 
of inhomogeneity. 



IV. RESULTS AT FINITE T 

The SC transition temperature Tc of a lattice with 
an inhomogencous pattern of interacting sites can also 
be larger compared to a uniform interaction distribution 
on the same lattice. In this section, we investigate the 
conditions under which inhomogeneity in any form can 
lead to the enhancement of Tc as a function of |C/| and 
n at different concentration values /. In Fig. [7l panel 
(a) presents the variation of the average order parame- 
ter Ainhom and Auniform as functions of \U\ for / = 0.25, 
n = 1.35 < 1.5 = n*, and T = 0. We pick the value 
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FIG. 6: (Color online) The magnitude of the averaged or- 
der parameter A for the inhomogeneous patterns along the 
enhancement boundary (r = 1) in Fig. [T] as a function of 
the average electron occupation n for / = 0.25 (panel (a)), 
/ = 0.5 (panel (b)) and / = 0.75 (panel (c)). 



of \U\ = 6t, the largest in panel (a), at T = and plot 
both Ainhom(r) and Auniform(T) as functions of T in the 
inset inside panel (a). As seen in panel (a), for \U\ — 6t, 
Auniform(r = 0) has already exceeded all its inhomo- 
geneous counterparts appreciably. Nevertheless, in the 
inset inside the same panel, Tc for the inhomogeneous 
patterns are still larger than their uniform pattern coun- 
terpart at \U\ ^ 6t indicating the strong enhancement of 

For the uniform pattern at all values of \U\ and n, 
we find Tc in good agreement with the BCS prediction. 
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FIG. 7: (Color online) Panel (a): (main) The variation of the 
averaged order parameter A for the uniform and three dif- 
ferent inhomogeneous patterns of checkerboard, stripes and 
random at / = 0.25, n = 1.35 < n* and T = as a function 
of 1 17 1 . The inset presents the evolution of these order param- 
eters as a function of temperature for \U\ = 6t (The largest 
at T = 0). Panel (b): (main) The same results as in panel (a) 
for n = 1.5 = n* and T = 0. The inset shows the evolution of 
order parameters against temperature for two different values 
of \U\ = 2t (solid line with filled symbols) and 41 (dashed line 
with open symbols) as indicated in the T = results by the 
dashed lines. 



kbTc « (A(0)|[/|)/1.76, as expected from our mean field 
approach. In panel (b) in Fig. [7|for n = 1.5 ~ n*, how- 
ever, for \U\ = 2t for which Auniform(r = 0) is slightly 
larger than Ainhom(T' — 0), we find that Tc for all in- 
homogeneous patterns (except the checkerboard whose 
AinhomCT" = 0) = at |i7| = 2t) are also slightly larger 
than the uniform pattern Tc as shown in the inset of 
the same panel. At \U\ = At for which Auniform(T = 0) 
becomes noticeably larger than Ainiiom(r = 0), as pre- 
sented in the inset inside panel (b), Tc for the uniform 
pattern also becomes larger than its inhomogeneous pat- 
tern counterparts. In other words, the enhancement of 
Tc is rather weak when n > n* compared to n < n* 
and ceases to persist as |f7| increases. In panel (b), the 
checkerboard pattern has a vanishing average order pa- 
rameter at both T = and finite T. For stripes, Ajnhom 
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FIG. 8: (Color online) Panels (a), (b) and (c) refer to / = 0.5 
and n — 0.9 < n* , n = 1 = n* and n = 1.2 > n* respectively. 
Two different values of |?7| have been chosen in panels (b) and 
(c) at finite temperature as shown by the dashed lines in the 
T — results. In panel (b), for better visibility, results for 
stripes at \U\ = 3.5t (dashed line with open diamonds) are 
shown in magenta inside the inset. 
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FIG. 9: (Color online) Panels (a), (b) and (c) refer to f — 
0.75 and n = 0.4 < n* , n = 0.5 = n* and n = 0.65 > n* 
respectively. Two different values of |C/| have been chosen in 
all panels as shown by the dashed lines in the T = results. 
In panel (c), results for checkerboard at \U\ = 1.5t (dashed 
line with open diamonds) are shown in cyan inside the inset 
for better visibility. 



starts with a finite value at T = and gradually vanishes 
as T increases. For the random pattern, Ainhom starts at 
a value very close to zero at T = 0. However, at \U\ = At 
as T increases, there is a slight rise in the magnitude 
of Ainhom over some finite temperature window before it 
completely vanishes at high enough T. 



In Fig. [51 panel (a) presents similar results for / = 0.5 
at n = 0.9 < 71* and \U\ — At. As illustrated in the 
inset of the same figure, even at \U\ = At for which 
Ainhom = 0) < Auniform(T = 0), aU inhomogeucous 
patterns lead to larger Tc compared to the uniform pat- 
tern. In panel (b), for n = 1 = n*, at \U\ = 2.At, 
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FIG. 10: (Color online) Anomalous behavior of the averaged 
order parameter A as a function of temperature T at / = 0.75 
and n = 0.5 as presented in Fig[9] panel (b) (inset), for stripes. 
Panel (a) illustrates the evolution of the DOS as a function of 
temperature for \U\ = 1.86t. Panel (b) demonstrates how the 
local Ai on any of the individual four sites inside the block 
shown in Fig[T] panel (b) vary as a function of temperature 
T. Panel (c): The evolution of the local occupation number 
m on any of the individual four sites inside the block shown 
in Fig[l] panel (b) as a function of temperature T. 



AinhomlT" = 0) < Au„ifotm(7' = 0). However, except for 
the checkerboard, we still find an increase in Tc due to in- 
homogeneity. Similar to Fig. [7]panel (b), for both striped 
and random patterns, we also find a gradual increase in 
Ai„hom(T') as T increases before it totally vanishes at 
high enough T. The enhancement of Tc continues to per- 



FIG. 11: (Color online) Same as Fig. 10, but for the random 
pattern. 



sist up to \U\ = 3.5t, although Ainhom(r = 0) « for all 
inhomogeneous patterns, due to the gradual increase of 
Ai„hom(r) with temperature. Panel (c) corresponds to 
n = 1.2 > n* and it can be observed that for \U\ = 1.5t, 
for stripes only, Tc is increased. However, at slightly 
larger |C7| = 1.8t, Tc for the uniform pattern significantly 
dominates the inhomogeneous pattern ones and no grad- 
ual increase in Ainhom(?^) occurs as T increases. 

Fig. EJa) has / = 0.75 and n = OA < n*. Similar to 
the behavior observed for / = 0.25 and 0.5 when n < 
n*, up to the largest value of \U\ — 3t, inhomogeneity 
significantly increases Tc as illustrated in the inset of the 
same panel. This occurs despite the fact that Ainhoin(T = 
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0) < Au„iform_(T = 0). For n = 0.5 = n*, even at |i7| = 
1.86i where Aiiihom(2^ = 0) « 0, for both the striped 
and random patterns we again find an increase in Tc due 
to inhomogeneity. In this case, the gradual increase in 
(T) as a function of T is further pronounced until 
Tc totally vanishes for these two patterns at high enough 
T. Panel (c) corresponds to n = 0.65 > n* and similar 
to / = 0.25 and / = 0.5 cases when n > n* , for large 
enough |i7| = 1.5t (slightly larger than |i7| = 1.25t in the 
inset) inhomogeneity no longer yields larger compared 
to the uniform. (See inset.) 

The anomalous increase of Ainhom as a function of T 
at n = n* for / > 0.5 is an actual feature and is believed 
to be related to the gradual destruction of the charge 
ordered phase due to temperature leading to an inter- 
mediate SC phase. In Fig. [TO] corresponding to Fig. [SI 
panel (b) for n — 0.5, / — 0.75 and \U\ = 1.86t, panel 
(a) presents the DOS results at several T values within 
the temperature window of the inset in Fig. [9| panel (b). 
At r = 0, the gap in the DOS is barely non-zero at 
a; = (Fermi energy) indicating a charge ordered phase. 
By increasing T to 0.6t the gap widens towards super- 
conductivity consistent with the behavior shown inside 
the inset in Fig. [9l panel (b). By further increasing T, 
the gap begins to shrink due to the destruction of super- 
conductivity by temperature until it entirely vanishes at 
T = 1.65t in agreement with the results inside the inset 
in Fig. [9l panel (b). 

Panel (b) in Fig. [TOlillustrates the variation of the local 
order parameter on all four different types of sites for 
stripes at / = 0.75. Only Ai along the line of interacting 
sites varies appreciably as T increases. Panel (c) of Fig. 
1101 shows the local occupation numbers on all the four 
different types of sites for stripes at / = 0.75 (Fig. [TJ 
panel (b)). By increasing T, while the occupation of the 
interacting sites gradually drops from n = 2 to n ~ 1, an 
indication of the destruction of the charge ordered phase, 
the occupation of the non-interacting sites grows leading 
to additional mobility and overlap of the Cooper pairs 
and therefore enhancing Tc- 

Fig. [TI] panel (a) shows the DOS results for the ran- 
dom pattern for n = 0.5, / = 0.75 and \U\ = 1.86t 
depicted in Fig. [9l panel (b). Similar to stripes, over the 
window of temperature studied in panel (a), the gap in 
the DOS initially grows up to T « t and then gradually 
drops to zero slightly beyond T = l.Qt to indicate the for- 
mation of a metallic phase as the temperature destroys 
superconductivity. Panel (b) in Fig. [TT]also shows signif- 
icant variation of the local order parameter only on the 
interacting sites on average. Hence, similar to stripes, 
temperature acts against the more localized charge or- 
dered phase allowing for the Cooper pairs to move and 
overlap more freely and consequently the enhancement 
of Tc- Panel (c) also confirms that the charge ordered 
phase is indeed destroyed by temperature allowing the 
non-interacting sites to become more occupied resulting 
in an intermediate SC phase. Metallic behavior is estab- 
lished at sufficiently high temperature. 



It is also noteworthy that this anomalous increase in 
Ainhom was not observed for the checkerboard pattern for 
which there is no superconductivity at any temperature 
when n = n* . One possible explanation is that due to 
the particular geometry for the checkerboard, localized 
Cooper pairs are further apart from one another com- 
pared to the striped and random patterns. As mentioned 
earlier, when n = n* we have two electrons per interact- 
ing site. For the checkerboard lattice, this leads to a 
pair localized in the interacting sites with empty nearest 
neighboring sites. Therefore, the effect of temperature is 
minor in causing further overlap among the pairs before 
it totally destroys them (especially at / = 0.5 and 0.75 
as shown in Fig. [TJ panel(a)). 

A similar increase in Tc upon introducing a checker- 
board pattern is found in the MCMF calculations as well, 
arising from the loss of long-range phase coherence. This 
is particularly significant because the MCMF incorpo- 
rates the subtle nature of the SC transition in 2D dis- 
cussed earlier. We have also independently confirmed 
that our conclusions and arguments equally apply for a 
model with nearest-neighbor attraction, leading to a d- 
wave SC close to half-filling, which refiects the cuprates' 
phenomenology more truthfully42.i2i 

We have found that upon introducing inhomogeneity 
into the pattern of interacting sites on a lattice. Tc can 
be increased over a wide range of \U\ as long as n < 
n* even if Ainhom (T = 0) < Au„ifo„„(r = 0). When 
n > n* , however, for sufficiently large \U\, inhomogeneity 
acts against superconductivity and therefore Tc becomes 
smaller compared to the uniform pattern of interacting 
sites. The case oi n — n* is anomalous. The charge 
ordered phase established at large enough |i7| values at 
T — gradually vanishes as T increases. Ainhom first 
increases with increasing T, leading to a SC phase (at 
least for large enough / values and lower symmetrical 
inhomogeneity patterns), and then vanishes, indicating a 
metallic state. 



V. SUMMARY AND DISCUSSION 

In summary, we have shown that for the attractive 
Hubbard model on a square lattice, there is a signifi- 
cant range of electron doping and interaction strength 
over which the average superconducting order parame- 
ter is larger for a lattice with inhomogeneous patterns of 
interacting sites than a uniform distribution of these in- 
teracting sites at a constant interaction strength per site. 
We have presented the phase diagrams for three different 
inhomogeneous patterns: checkerboard, stripes and ran- 
dom and also three different values for the concentration 
of the non-interacting sites. Apart from a few particu- 
lar features, the overall physics illustrated in the phase 
diagrams is pattern independent. As we vary the mean 
interaction strength \U\ and the doping level n at zero 
temperature, we have verified the existence of at least 
three different phases, namely (i) superconducting, (ii) 
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insulating due to the charge order phases, and (iii) metal- 
Hc. Our findings and claims are strongly supported by 
studying the behavior of a variety of quantities computed 
in this work all consistently corroborating one another. 

The enhancement of the average order parameter for 
the inhomogeneous interacting site patterns is due to the 
proximity effect, i.e., the tunneling effect of the Cooper 
pairs from the interacting sites leading to finite order 
parameter values on neighboring sites. This conclusion 
is supported by the effect occurring at weak coupling, 
where the coherence length is large, rather than in the 
strong coupling regime of preformed pairs. Agreement 
between the BdG results and the MCMF calculations 
justifies the application and conclusions of the BdG ap- 
proach within the small \U\ regime. Our calculations 
also clearly confirm that an inhomogeneous interaction 
potential can lead to the increase in the phase transition 
temperature Tc over a wide range of n, and \U\ for vari- 
ous / values. Counterintuitively, as long as n is less than 
or equal to twice the fraction of interacting sites, this in- 
crease in Tc continues even for values of for which the 
order parameter is larger for the uniform pattern than for 
inhomogeneous patterns at T = 0. 

One possible explanation takes into account that in 
this weak coupling parameter regime, Tc is a supralin- 
early increasing function of U. In such a case, it may be 
that in the inhomogeneous system the sites with larger 
U produce a nonlinear enhancement relative to Tc of the 
uniform system and, through the proximity effect, drag 
the non-interacting sites along with them. This trend 
changes when n exceeds twice the number of interacting 
sites (i.e. some electrons must occupy non-interacting 
sites), for which at large enough \U\ values inhomogenc- 
ity fails to increase Tc over that of the uniform pattern. 
The n = 2(1 — f) = n* case for sufficiently large / val- 
ues and less symmetric inhomogeneous patterns (such as 
stripes and random as opposed to the checkerboard) is 
anomalous as it shows the enhancement of Ajuhom (T) as 
temperature increases. 

It is even more surprising to find that a system which is 
non-superconducting (charge ordered insulator) at T = 
can become superconducting upon increasing T for a 
finite window of temperature before turning metallic. 
This anomalous behavior was shown to be related to 
a crossover from a charge ordered insulating phase for 
n = n* at large enough \U\ values to an intermediate 
SC phase upon increasing T before entering the metallic 
phase at sufficiently large T. 

We wish to emphasize that in this article we have fo- 
cused on the enhancement (or not) of the pairing ampli- 
tude (our A defined in section|TT|, rather than its product 
with the local interaction strength —\Ui\ which is more 
directly related to the gap but contains less information, 
and gives less insight, because it automatically vanishes 
on any site without interaction. Thermodynamic mea- 
surements would probe quantities which include the en- 
ergy scale, such as the specific heat or superfiuid rigid- 
ity which our results may not have direct implications 



to. However, for an inhomogeneous system being a mix- 
ture of different phases, defining an average SC gap is 
not trivial. Thus, the thermodynamic properties of in- 
homogeneous superconductors will not necessary exhibit 
the same behavior as their homogeneous counterparts. 
It has been shown that the rise of the specific heat in 
inhomogeneous superconductors obeys a power law be- 
havior as opposed to exponential in homogeneous ones 
using the attractive Hubbard model with random inter- 
acting sitesj^ Also, the superfiuid density and stiffness in 
general decrease due to the presence of disorder (22, Never- 
theless, lower superfiuid density does not necessarily lead 
to lower Tc as according to the Anderson theorem, a non 
magnetic impurity should not affect the Tc and therefore 
thermodynamic properties of a s-wave superconductor. 
We have presented clear evidence for the enhancement 
of Tc which does have a direct experimental implication. 
Thermodynamic properties of inhomogeneous supercon- 
ductors are very rich in physics and a great deal of con- 
tributions and new ideas in this area are yet to appear. 

While the attractive Hubbard Hamiltonian obviously 
does not incorporate many of the features of high Tc su- 
perconductors (notably the symmetry of the pairing), the 
model has been shown to provide useful insight into some 
of their phenomenology, for example the spin-gap?^ It is 
therefore tempting to speculate that our results concern- 
ing inhomogeneity may have similar connections. Specif- 
ically, earlier ARPES data^ suggests that the under- 
doped phase of LSCO (La2-a;Sra;Cu04) consists of SC 
clusters, embedded in the AF host. In such a system, in- 
homogeneous gaps appear naturally and our results here 
indicate that the SC transition is in fact determined by 
the largest gap values rather than the much smaller gaps 
found at phase boundaries, as one might naively think. 
This renders the SC phase more stable than it would 
otherwise be, and also simplifies the description of these 
systems. 

It is worth emphasizing that in most situations, in- 
homogeneities reduce values of order parameters^ and 
critical temperatures, even when comparisons are made, 
as they are in this article, to homogeneous systems with 
the same average value of all parameters. This is true, 
for example, for classical site diluted Ising models, where 
the ferromagnetic J is increased to compensate for ab- 
sent sites, and quantum models like the boson Hubbard 
model where random chemical potentials monotonically 
decrease and ultimately destroy superfiuidity.^ -'S2 An ex- 
ception is the increase of TnccI by randomness reported 
in DMFT studies of the repulsive model®^ and recently, 
the SC Tc in XY model Hamiltonian with certain types 
of inhomogeneous patterns for the coupling constant be- 
tween spins sitting on nearest neighboring sitesi^. 
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